Aquatic Conservation: Marine and Freshwater Ecosystems

The Peruvian anchovy (Engraulis ringens) is a small pelagic fish. It is the largest singlespecies fishery in the world and it is used mainly for the production of fishmeal and fish oil. The Peruvian hake (Merluccius gayi) is a predatory demersal fish that is valued for direct human consumption, mainly as frozen products. In this study, I compared the response of both species to certain fishing mortality scenarios using a climatological set-up of the multispecies model OSMOSE for the northern Humboldt Current System. I observed that hake landings benefit from a decreased fishing pressure. In addition, the resilience of the anchovy fishery may increase if a higher threshold for their minimum catch size was implemented. The results of this study provide insights into better management strategies that could be beneficial for the two species based on their life-strategies.

Environmental concerns and the need to mitigate the effects of climate change have driven the adoption of renewable energy sources, such as offshore wind farms (OWFs), making it essential to assess their overlap with seabird habitats to estimate potential impacts. This study is the first to estimate the vulnerability of coastal and seabirds to OWFs at a national scale by integrating biological trait data from the literature, spatial occurrence records from citizen science platforms, and expert knowledge (elicited from 23 specialists) into established quantitative risk frameworks. Collision risk and susceptibility scores were calculated at the family and subfamily levels, as well as disturbance and displacement indices by species and susceptibility by state. The groups with the highest collision risk and susceptibility were Fregatidae (0.354; 0.486), Sulidae (0.199; 0.419), Sterninae (0.190; 0.384) and Diomedeidae (0.190; 0.416). The extreme risk for Fregatidae is driven by their specific flight characteristics, with impacts concentrated in the highly frequented southeast and northeast regions. Regarding disturbance and displacement, Limnodromus griseus exhibited the highest national score (MMA: 0.098), while Diomedea dabbenena had the highest global score (IUCN: 0.092), highlighting the critical need to integrate both regional and global conservation statuses. Spatially, the state exhibiting the highest overall vulnerability to OWF impacts was Rio Grande do Sul (23.78; 10.69), a critical hotspot characterized by high species richness and oceanographic productivity. These findings establish an indispensable national baseline for Marine Spatial Planning, reinforcing the need for targeted pre-installation assessments, the identification of sentinel species, and spatially differentiated management during environmental licensing processes.

Sea turtle conservation relies largely on the identification of ecologically significant marine areas. Green sea turtles (Chelonia mydas) are resident in Aotearoa New Zealand waters, yet much of their spatial ecology remains undefined. To address this, satellite telemetry was used to map juvenile and sub-adult green sea turtle spatial movements, residency patterns, and habitat associations. Our findings provide the first evidence of individual fine-scale habitat use and long-term residency in New Zealand waters. Data will contribute to regional and global efforts to identify ecologically significant sea turtle habitat, including Important Marine Turtle Areas (IMTAs) and Blue Corridors for Turtles (BC4T), and help advance conservation efforts.

The decline in fishery yields in Japan's enclosed coastal seas has been accompanied by a reduction in phytoplankton biomass, suggesting a widespread decrease in ecosystem productivity. Oligotrophication in Mikawa Bay has caused substantial declines in phytoplankton biomass and Manila clam (Ruditapes philippinarum) production. To assess whether external nutrient enhancement could stimulate phytoplankton growth, a controlled demonstration experiment using treated effluent from two sewage treatment plants (WWTPs) was conducted (spanning fiscal years 2022 to 2024) based on analysis using a three-dimensional non-hydrostatic model. Although effluent DIP and nitrogen concentrations were increased, model results showed that WWTP-derived nutrients represent only a minor portion of the total nutrient budget. Field observations revealed limited increases in surface water phosphorus, strong seasonal variability, and only temporary rises in chlorophyll-a and clam harvests. Monthly TP data demonstrated that seasonal DIP pulses, especially spring phosphorus bursts, dominate nutrient dynamics and overshadow external loading effects. Laboratory experiments further showed that phosphate is efficiently fixed by Fe(III) under oxidizing conditions and rapidly released during reductive dissolution of [FeOOH-HPO4]2- complexes, confirming iron as the principal regulator of internal P cycling. Because current numerical models omit iron-driven internal cycling, they cannot accurately predict bioavailable P or ecosystem responses. Overall, the study indicates that localized nutrient enhancement is insufficient to counter long-term oligotrophication, and effective management must integrate nitrogen limitation, internal Fe and P cycling, oceanic exchange, and seasonal variability in Ise and Mikawa Bays.

The identification and use of stress-tolerant coral in outplanting efforts has been proposed as a method of restoring and bolstering reef health under increasing environmental change. Prior research has focused on understanding survival in cross-reef transplantation, but evidence of long term resiliency of outplanted corals remains limited. Here, we report data at a five-year timepoint of a reciprocal transplant study of the stress-tolerant Caribbean corals Pseudodiploria strigosa and Siderastrea siderea between a nearshore and offshore site in Belize. Extreme marine heat wave conditions in August of 2024 revealed species and site-specific bleaching patterns. Through assessments of mortality, energy reserves, and chlorophyll concentrations, P. strigosa appear to be more limited in long-term acclimatization potential to the nearshore environment while S. siderea maintain acclimatization capacity displayed 17 months post transplant to nearshore. This work conveys information relevant to managers in selecting species for restoration efforts based on short-term growth and long-term resilience in native versus non-native sites, offers insight into potential mechanisms driving these trends, and takeaways on site specific management decisions for successful coral outplanting success.

1. Protecting and restoring corals is critical to conserve marine biodiversity and ecosystem services in the Anthropocene. Knowledge of coral larval dispersal dynamics can be used to identify important sites for coral metapopulation connectivity and to spatially optimise coral conservation efforts. However, we know little regarding the larval dispersal of many coral species and our understanding of the relationships between coral larval dispersal and life-history traits (e.g., reproductive mode, colony morphology), phylogenetic relationships, and biogeography is incomplete. These knowledge gaps impede coral conservation. 2. In this research, we systematically reviewed the species, regions, and methodologies of studies of coral larval dispersal. Additionally, we conducted a meta-analysis to quantify the heterogeneity in reported coral larval duration values and identified the roles of life-history traits in driving that heterogeneity. 3. We found that various methods (experimental, biophysical modelling, and population genetics) were used to study 221 coral species across many regions. Soft and temperate corals received comparatively lower research attention. Our meta-analysis reveals that reproductive mode and colony morphology explain significant heterogeneity in larval settlement times, and that closely-related species have more similar larval settlement times. 4. Our findings suggest that reproductive mode and colony morphology could be used to predict larval duration in under-studied coral species, and that larval duration could be approximated from closely-related corals when trait data is unavailable. Our results support strategies focused on modelling differing larval durations across corals with different life-history traits to predict realised dispersal in under-studied coral species and identify conservation priority sites.

Introduction: Despite decades of investment, stream restoration often fails to recover diverse benthic macroinvertebrate communities, as restored reaches tend to remain dominated by pollution-tolerant taxa. This shortfall in recovery is commonly attributed to either poor habitat quality or dispersal limitation of aquatic macroinvertebrate taxa, yet these hypotheses are rarely tested directly. Objectives: We conducted a manipulative field experiment in the eastern Piedmont region (Maryland, USA) within the Chesapeake Bay Watershed to evaluate whether restored streams can support pollution-sensitive taxa when dispersal barriers are removed. Methods: We translocated whole macroinvertebrate communities from three nearby reference streams into three restored streams using standardized habitat substrates and quantified macroinvertebrate survival after 28 days. We also assessed how effective leaf vs rock substrates were for accumulating diverse and abundant macroinvertebrate communities for translocation. Results: Across all sites, 77% of translocated individuals and 16 of 22 taxa sensitive and moderately sensitive to pollution persisted under restored stream conditions, demonstrating that restored reaches can support short-term survival of pollution-sensitive taxa that typically fail to recover. Substrate comparisons revealed that leaf packs accumulated higher macroinvertebrate abundances, whereas both substrates accumulated macroinvertebrate communities with similar taxa richness, EPT richness, and Shannon Diversity. Conclusions: This work suggests that dispersal limitation, rather than unsuitable habitat, may be an underappreciated constraint limiting macroinvertebrate recovery in restored streams. These findings challenge the assumption that habitat restoration alone is sufficient for biological recovery and highlight the potential for whole-community translocations to accelerate recovery and ecological uplift as part of stream restoration practices. Implication for Practice: Our study highlights the need to consider habitat connectivity and dispersal pathways during stream restoration designs. Restoring streams to meet habitat requirements of pollution-sensitive macroinvertebrate taxa may not be sufficient for their populations to recover if these taxa are unable to reach the restored stream. In scenarios where restoring connectivity to source populations is infeasible, whole community translocations may offer a novel restoration tool to overcome dispersal barriers and accelerate biological recovery in restored streams.

Diverse life history strategies allow species to spread risk across a mosaic of habitat conditions. In Chinook Salmon (Oncorhynchus tshawytscha), particularly in the Central Valley of California, this is expressed through varied adult run timings and juvenile outmigration strategies. Anthropogenic impacts, particularly dams, disrupt these fundamental life stage transitions, often forcing managed populations to adapt to homogenized, less stochastic hydrologic regimes. This is exemplified in the Sacramento River, where the conservation of endangered Sacramento River winter run Chinook salmon requires balancing complex water operations with the ecological needs of a population confined to a short stretch of suitable habitat downstream of Shasta Reservoir. Using a 23 year dataset (2002 to 2024), we evaluated the combined effects of flow, temperature, and spawner abundance on juvenile production and life history expression in this habitat. We found that flow and spawner abundance best predicted juvenile abundance at Red Bluff Diversion Dam while temperature had considerably less support. The proportion of juveniles migrating as smolts exhibited significant density dependence that was negatively correlated to both female spawner abundance and peak flows. These results suggest focusing on temperature management alone may be insufficient. Effective management strategies should consider flow variability and habitat restoration to facilitate varying migration strategies and expand upstream rearing capacity. By addressing these physical and hydrologic constraints, managers can better support the full suite of life history strategies necessary for the resilience of winter run Chinook salmon.

Reef fish assemblages in the Veracruz Reef System National Park (VRSNP) are strongly structured along a pronounced north-south gradient, where environmental conditions, urban proximity, and protection levels interact to shape functional diversity and ecosystem functioning. This gradient creates contrasting ecological conditions between northern and southern reefs, providing a natural framework for evaluating how spatial variation influences resilience. Using a stationary visual census across ten reefs distributed latitudinally, we quantified functional entities based on six ecological traits and calculated functional richness, redundancy, over-redundancy, vulnerability, specialization, and originality. A two-way permutational ANOVA and redundancy analysis were used to detect spatial and temporal patterns and to link functional indices with environmental features. The southern sector displayed the most intense functional shifts. Nearshore southern reefs, strongly exposed to sedimentation, nutrient loading, and physical disturbance, exhibited reduced functional richness, fewer functional entities, and high functional vulnerability, indicating strong environmental filtering. Over time, these reefs shifted toward generalist-dominated assemblages with rising redundancy, suggesting functional convergence. In contrast, southern offshore reefs supported higher taxonomic and functional richness, buffered by greater habitat complexity and distance from shore, although they showed declining redundancy, likely reflecting sensitivity to regional climatic disturbances. In the northern sector, the gradient was moderated by effective protection, which reshaped functional configurations. Nearshore protected reefs such as Isla Sacrificios and Isla de Enmedio maintained high redundancy, low vulnerability, and high biomass of large-bodied fishes, resembling offshore systems rather than disturbed southern nearshore reefs. The northern offshore reef (Isla Verde) shared some of these attributes but exhibited temporal declines in richness, signaling emerging vulnerabilities along the gradient. Functional originality and specialization aligned with benthic states across the gradient, particularly where macroalgal or calcareous substrates prevailed. Overall, the north-south gradient is the primary axis of functional organization in VRSNP, revealing that resilience emerges from the interplay among reef geomorphology, urban influence, and governance effectiveness.

Okanagan sockeye salmon have shown a strong recovery over the past three decades, but the contribution of environmental conditions and management actions to this recovery has not been determined. We developed a hierarchical Bayesian model linking spawner abundances to fry emergence observations to evaluate how flow and water temperature during adult upstream migration and incubation periods influence spawner-to-fry production. Production declined by 62% in years with earlier onset and longer duration of periods when warm river temperatures created barriers to adult migration, indicating that migration conditions can have carry-over effects on reproductive success even when spawners successfully reach the spawning grounds. Additionally, production was 46% lower when flows fell above or below the management target range designed to minimize scour and desiccation during egg incubation. In contrast, we found little support for effects of pre-spawning flows or incubation temperature on production. These findings provide quantitative evidence that current flow management is increasing freshwater spawner-to-fry production and has likely contributed to the Okanagan sockeye recovery. Effective flow management may therefore help offset climate-driven impacts that increasingly threaten this population.

Coral skeletal integrity underpins reef habitat complexity and coastal protection, yet it remains unclear whether corals persisting in chemically variable marginal habitats, such as mangrove-associated reefs, maintain sufficient mechanical resistance to withstand physical disturbance. We combined micro-computed tomography (micro-CT) morphometric analysis with compression testing to quantify skeletal microarchitecture and compressive mechanical performance of Pocillopora sp. skeletons collected across a habitat gradient at Low Isles, Great Barrier Reef, from mangrove-associated sites to exposed reef sites. Baseline skeletal properties varied among sites, with the Proximal Reef site having the densest and least porous skeletons, while mangrove-associated sites had the most porous skeletons. Similarly, Proximal Reef skeletons exhibited two-to-three-fold greater compressive stiffness and peak force than mangrove-associated skeletons. Percentage object volume was a stronger predictor of compressive stiffness than surface area-to-volume (S/V) ratio, indicating that bulk skeletal density better captured load-bearing capacity than surface architectural complexity. Structural failure mode also differed markedly among habitats: mangrove-associated skeletons underwent catastrophic fragmentation of their open, porous framework, whereas Distal and Proximal Reef skeletons failed through discrete crack-plane propagation while largely maintaining overall structural coherence. These results indicate that the physiological stress tolerance associated with mangrove habitats occurs alongside reduced mechanical resistance and more destructive failure modes, suggesting a potential trade-off between chemical stress tolerance and structural investment. Assessments of resilience in marginal coral populations should therefore consider skeletal mechanical performance alongside physiological tolerance, particularly as climate change intensifies storms and alters hydrodynamic disturbance regimes.

Hydropower development has emerged as a major driver of global environmental change, causing significant ecological damage to freshwater ecosystems and threatening aquatic biodiversity. However, the impact of cascade dams on fish communities in large river systems remains poorly understood, particularly in biodiversity hotspots such as China. Here, we used environmental DNA (eDNA) metabarcoding to assess the spatiotemporal patterns of fish diversity and community structure in the Ganjiang River, the largest tributary of the Yangtze River, in relation to cascade hydropower development. We collected water samples from 19 sites along the main stem and major tributaries in both wet and dry seasons, and analyzed the eDNA data using high-throughput sequencing and advanced bioinformatics pipelines. Our results revealed significant seasonal and spatial variations in fish diversity and community composition, with lower species richness and distinct assemblages in the impounded sections upstream of the dams compared to the downstream reaches. The fish communities were strongly influenced by environmental factors such as water temperature, pH, dissolved oxygen, and nutrients, which exhibited complex interactions with dam operation and land use practices. Spatial scale analyses further indicated that the impact of dams on fish diversity was most pronounced at the local scale, but also detectable at the regional scale. Our findings provide novel insights into the ecological consequences of cascade hydro-power development on fish communities and highlight the importance of considering spatial scales and environmental heterogeneity in assessing the impacts of human activities on freshwater biodiversity. Our study demonstrates the power of eDNA metabarcoding as a cost-effective and non-invasive tool for biodiversity assessment and conservation planning in large river systems facing multiple stressors and knowledge gaps.

Management of marine protected areas (MPAs) require long-term monitoring and evaluation that is often limited by the cost and effort of conventional visual based surveys. We evaluated a community science-based eDNA approach to MPA monitoring over two years in the Point Dume State Marine Reserve, the Point Dume State Marine Conservation Area, and two paired unprotected reference sites. Across the sampling period, we detected 119 species of fish and marine mammals, including 19 of the 22 regional management-priority species for Southern CA MPAs. Protected and unprotected comparison sites were broadly similar in overall assemblage with 93 species shared across all sites. Multivariate analyses revealed site-level differences, although season and year explained more variation than protection status. Species presence also shifted during a La-Nina transition, highlighting the importance of long-term monitoring under different climatological contexts. Community scientists greatly expanded the temporal and spatial scope of eDNA sampling and generated high-quality data. Subsampling analysis showed that quarterly sampling, with three sample replicates and sequencing depths between 700k-990k, recovered 95% of observed biodiversity of the full time series. These results demonstrate that a community science-based eDNA program could provide a standardized, scalable, low-cost, and management-relevant approach for long-term monitoring across California's MPA network.

Understanding the consequences of global change for our ecosystems requires reliable, accurate, and scalable tools for detecting ecological change. To this end, recent efforts have focused on scaling up the quantification of demographic processes from individual organisms (the fundamental ecological unit) to the larger spatial scales they inhabit to yield more precise and informative estimations of ecological change, especially in response to disturbance. Here, we describe a transferable workflow that combines robust underwater photogrammetry techniques with machine-learning-assisted image analysis to quantify demographic changes in coral colonies (growth, stasis, partial mortality, and complete mortality) through time on the reef landscape. We then demonstrate several analytical applications of this workflow that can be broadly applied to sessile and/or clonal organisms in other ecosystems. First, we tracked the fates of 3,162 coral colonies through a marine heatwave that caused substantial mortality. Second, our approach enabled detection of size- and taxon-dependent patterns in demographic changes from before to after the heatwave. Finally, for a subset (425) of these colonies, we further demonstrated how our method can be used to track fates across multiple years, including a year of change before and after the heatwave took place, revealing nuanced patterns of survival and mortality, as well as changes in size structure that occurred before and after the disturbance. Our workflow offers an important advance for increasing the accuracy of change detection on coral reefs and other ecosystems by linking demographic changes in individual organisms with spatially extensive image-based surveys through time.

The genus Phakellia comprises predominantly fan-shaped deep-sea sponges, but its species diversity and boundaries remain unresolved. Here, we use an integrative approach combining morphological and molecular data to investigate the Phakellia robusta species complex and describe two new species, Phakellia pallida sp. nov. and Phakellia pusilla sp. nov. Overlap in external and spicular morphology makes these taxa challenging to identify. Molecular analyses based on mitochondrial (cox1) and nuclear (28S) markers support the recognition of the three taxa as distinct species, and reveal their close evolutionary relationship. The two new species co-occur in the NE Atlantic Ocean and the Mediterranean Sea, with overlapping depth ranges, although P. pallida sp. nov. is generally found at greater depths. Our results also suggest that several previous records of P. robusta from the Mediterranean and southern parts of the NE Atlantic may instead correspond to either P. pallida sp. nov. or P. pusilla sp. nov., indicating that its distribution requires reassessment. In addition, a group I intron previously reported in the cox1 gene of P. robusta was polymorphically detected at the same insertion site in some, but not all, individuals of P. pusilla sp. nov., but it was absent from P. pallida sp. nov. individuals. Overall, our results reveal pseudocryptic diversity within Phakellia and highlight the value of integrative approaches for resolving species boundaries and reassessing deep-sea sponge diversity and distributions.

This study assesses the spatial dynamics of heavy metals and associated trace elements in the water and sediment of the Slana River approximately one and half year after an environmental accident and river contamination caused by mining effluent. Using PCA and standardized environmental indices (CF, PLI), we modelled the distribution of contaminants across the macro-gradient of the river and within the micro-gradient surrounding the source of pollution. We found that although the situation in the water column has improved and stabilized, the sediment remains the primary reservoir of trace elements, including toxic arsenic, with the highest concentrations recorded near the source and in its immediate vicinity. Although the concentration of primary mining indicators decreases relatively rapidly along the stream, this is systematically offset by the accumulation of other elements, including sulphur and aluminium. As a result of these geochemical processes, the total pollution load index (PLI) remains statistically unchanged along the entire stream, but indicates a chronic, spatially uniform (PLI = 1.06-1.16) toxicological burden on the ecosystem.

Environmental nucleic acid (eRNA and eDNA) metabarcoding is a promising tool for biodiversity assessment as it supports low-cost, noninvasive biological surveys. A critical challenge for eRNA/eDNA analyses is accounting for nucleic acid contamination that can compromise the accuracy of biomonitoring results. The primary aim of the current study was to examine carcasses as potential sources of eRNA/eDNA contamination, a contribution that remains poorly understood. Using Japanese medaka (Oryzias latipes), we quantified the detectability and persistence of carcass-derived eRNA/eDNA in water through two complementary experiments: (i) a shedding study comparing eRNA/eDNA release from live and dead individuals and (ii) a decay experiment measuring temperature-dependent degradation from 4 to 34 {degrees}C. Dead fish released detectable amounts of both eDNA and eRNA, indicating that eRNA-based surveys cannot fully eliminate carcass-derived signals. Arrhenius-based analysis showed that eRNA degradation was more temperature-dependent than eDNA degradation. These results suggest that eRNA may improve temporal resolution for recent organismal presence under warm conditions, whereas carcass-derived eRNA may persist under cold conditions and complicate interpretation of positive detections. Overall, eRNA may help reduce carcass-related false positives, particularly in warm conditions, but should be paired with temperature-aware sampling strategies and temporal replication during colder periods.

Fledglings of burrow-nesting seabirds are vulnerable to disorientation by artificial light at night (ALAN) and grounding in coastal human settlements, where most individuals presumably perish if not rescued. Accurate estimates of the effects of stranding on seabird populations remain elusive due to several levels of uncertainty, such as the number of stranded seabirds estimated to remain undiscovered by rescuers. One species subjected to stranding is the Atlantic puffin (Fratercula arctica). We hypothesized that the detection probability of puffin fledglings by the Puffin Patrol rescue program, operating in Newfoundland, Canada, is affected by light level, concealment level, and search effort. We deployed realistic puffin decoys in locations of varying light and concealment levels throughout two sites that are searched nightly by rescuers and measured the amount of time and search effort required until each decoy was found. Approximately 78% more search effort was required to find a decoy in darkness (32.0 min) compared with illuminated conditions (18.0 min), and up to 100% more search effort was required to find a decoy placed in concealed conditions (35.5 min) than one in the open (17.5 min). After two hours of searching, on average over 10% of decoys were still not found, although each additional minute of search effort increased detection rate of an individual decoy by 0.6%. In combination with recent advances in our knowledge of seabird post-stranding behaviour, our study offers a first step in the development of correction factors which can be adapted to site- and species-specific requirements of other rescue programs around the world that rely on volunteers to search for stranded birds. Application of such correction factors will hopefully result in better estimation of mortality rates and their resultant effects on populations.

The transformative impacts of global bleaching events are now a ubiquitous feature of coral reefs, often resulting in taxonomically simplified coral assemblages that can persist long-term due to ongoing ocean warming, heatwaves, and localized anthropogenic stressors. Understanding to what extent these processes shape symbiont biodiversity continues to be a major knowledge gap. Here, we examined Acropora host and symbiont communities on Kiritimati (Republic of Kiribati), seven and a half years after the 20152016 El Nino-associated marine heatwave, which caused a near total loss of the genus in the atoll's shallow reefs. We quantified Acropora cover and growth-form composition at 37 sites spanning a gradient of local human disturbance and characterised Symbiodiniaceae communities using ITS2 DNA metabarcoding from 109 coral tissue samples at 13 of the sites. Corymbose and tabular Acropora growth forms remained evenly represented across the atoll and cover declined with increasing human disturbance. Symbiodiniaceae communities were dominated by Cladocopium with 93 % of colonies associated with a single ITS2 profile. Communities were homogenous, with low diversity and insignificant differences amongst growth forms, but slight differences across disturbance levels. DNA barcoding of the PaxC intron identified the majority of host colonies as most similar to Acropora granulosa (62 %) and Acropora hyacinthus (24 %). Although dominance by a single symbiont lineage is common among Acropora species across the Pacific, the low Symbiodiniaceae diversity observed at Kiritimati contrasts with the substantial within- and among-lineage diversity documented elsewhere. The depauperate Symbiodiniaceae communities observed in this study may thus represent a lasting legacy of the 2015-2016 disturbance, with the recovering Acropora assemblage representing an early stage of community reassembly drawing from a potentially depleted environmental symbiont pool.

The increased frequency and intensity of coral bleaching events expose reefs to disturbance-recovery cycles that threaten their long-term dynamics. Understanding the impact of thermal stress requires evaluating coral mortality and recovery rates after disturbances. By integrating two long-term monitoring datasets, we examined coral mortality and recovery patterns at 13 sites in Mo'orea, French Polynesia, following the 2019 bleaching event, the most severe heat stress on record for the island. The maximum accumulated heat stress (AHSmax) varied from 5.76 to 6.26 {degrees} C weeks and enabled a direct analysis of heat-stress-related mortality. One year after bleaching, Pocillopora mortality increased with AHSmax, reducing mean coral cover from 35% to 17%. Regardless of AHSmax, Acropora mortality was high at all sites (median: 84%), reducing mean cover from 5% to 1% and indicating greater heat tolerance of Pocillopora than Acropora. Within four years after the bleaching event, no recovery occurred, and coral cover further declined relative to pre-bleaching levels. This decline was stronger for sites exposed to more intense heat stress. The impact of colony size was assessed by tracking 3090 individual Pocillopora colonies over six years and estimating size-dependent mortality and growth with a Bayesian mixed model. Immediately after the bleaching event, large Pocillopora corals exhibited a higher probability of complete and partial mortality than smaller colonies, highlighting size-specific consequences of bleaching. Even three years after the heat stress ended, large Pocillopora corals remained stressed, as evidenced by the continuously high probability of partial mortality. The lack of recovery observed here may be related to increased spatial competition with macroalgae. The disappearance of large fecund colonies and high mortality of coral recruits could have catastrophic consequences for the resilience of Mo'orea's reefs.